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<N ■ Abstract 

h ; 

In the deconfinement phase of QCD quarks and gluons interact with the dense stochas- 
tic colour-magnetic vacuum. We consider the dynamics of quarks in this deconfinement 
■ phase using the Field Correlators Method and derive an effective nonperturbative inter- 

£NJ \ quark potential, in addition to the usual perturbative short-ranged interaction. We find 

the resulting angular-momentum-dependent interaction to be attractive enough to main- 
tain bound states and, for light quarks (and gluons), to cause emission of quark and gluon 
\£) , pairs. Possible consequences for the strong interacting quark-gluon plasma are briefly 

discussed. 
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£r ■ 1 Introduction 

O ; 

The picture of the strong interacting quark-gluon plasma (SQGP) seems to be adequate for 
explaining the recent data on ion-ion collisions pQ . It was found on the lattice [2] that colour- 
magnetic vacuum fields do not change across the phase transition, thus supporting the conjec- 
ture made in Ref. [3] that, above the temperature of deconfinement T c , the QCD vacuum loses 
its confining colour-electric part, while the colour-magnetic part remains intact. This idea can 
be economically expressed in the formalism of the Field Correlators Method (FCM) [I], where 
the Gaussian correlators of the colour-electric and colour-magnetic fields, ((Ei(x)Ej(y))) and 
((H i (x)Hj(y))) (((. . .)) denotes irreducible correlators), are parametrised through the correla- 
tion functions D E , D E , D H , and Df, respectively (see Ref. [I] for the details of the formalism^. 
Notice that among those only D E vanishes above the T c [21 E] ■ The usual (time-like) string 
tension o~e and the so-called spatial string tension a s are given by the double integrals from the 
corresponding correlators, 

°e = \J D E (Od% a s ^a H = l -J D H ($)d% (1) 

Both string tensions coincide below the T c but <te is expected to vanish in the deconfinement 
phase, whereas an remains nearly constant in the vicinity of the critical temperature, both 
below and above the T c jHlEj- Moreover, <jh grows at large T, as ^J(Jh oc Tg 2 (T) [6], signalising 
that D H grows as 0(T 2 g 4 (T)). At the same time it was conjectured in Ref. [7], and confirmed 
later on the lattice [21 E] , that the non-confining correlator Df does not vanish either above the 



1 The correlators D E and D H appear at the structure in these colour-electric and colour-magnetic cor- 
relators, respectively (see also Eq. ||5J)), whereas D 1 ' H parametrise the terms with derivatives. 
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T c and leads to strong interaction between quarks and gluons. In Refs. [H [9], bound states due 
to Df in quark and gluonic systems where found in qualitative agreement with the lattice data, 
Ref. [HI]. Many efforts based on the colour-electric type forces have been applied to clarify 
the dynamics of QCD above the deconfinement phase transition. Still, without the colour- 
magnetic forces, the dynamical picture of the SQGP is not complete since colour-electric forces 
cannot bind quarks and gluons for large angular momentum where colour-magnetic forces 
are most important. Thus the purpose of the present Letter is to clarify the matter and to 
study orbitally excited bound states of quarks above the T c . For the sake of simplicity, we 
consider only the confining correlators D E ' H and neglect the other two, D 1 ' . Furthermore, 
we study a strongly interacting system at the fixed temperature T, which can be either below 
or above the deconfinement temperature T c , and develop the standard Hamiltonian approach 
to it. We consider the situation when the colour-electric interaction is switched off above 
the deconfinement temperature and the dynamics of quarks and gluons is governed by the 
colour-magnetic forces only. Finally, we use the interaction under consideration to estimate 
the strength of the interaction in the SQGP by comparing the mean potential energy of light 
quarks (and gluons) in the plasma with their mean kinetic energy. We find the corresponding 
parameter to be of order <jh/T 2 . Since this parameter is large for quarks and even 9/2 times 
larger for gluons we conclude that SQGP is very strongly coupled and it should be viewed as 
a liquid, at least. 



2 QCD string and the spin— independent interaction in 
the quark— ant i quark system 

In this section we derive the effective spin-independent quark-antiquark interaction arising due 
the QCD string formation. 

Following the method proposed in Refs. [TTI [T2] we consider the gauge-invariant Green's 
function in the vacuum at temperature T: 

G(x u x 2 \y\, y 2 ) = (^outi^u x 2 )^L(yi, 2/2)), (2) 
where the wave functions of the initial and final colourless qq states are built with the help 
of the parallel transporter = Pexp i J* dz p A p (z) . Using the standard path integral 

approach and the Feynman-Schwinger representation for the single-quark propagator, one 
arrives at the Green's function ([2]) in the form [TTj : 

/OO I* poo f 

ds J (Dz)l yi e- K J drs J (Dz)l y2 e~ R (TrW(C)), (3) 

with the kinetic energies K = | J Q S drz^ and K = | J Q S drz 2 ^ and with the "winding path mea- 
sure" (Dz)™ y taking into account Matsubara periodic boundary conditions. All spin-dependent 
terms are neglected in Eq. ([3D — they will be restored below, in Section [31 The closed contour 
C runs over the quark trajectories and the dynamics of the system is defined by the averaged 
Wilson loop (TrW(C)). In the Gaussian approximation, one finds that [1] 

1 -(TrW(C)) = exp (-\ [ da^x) [ da Xp (x') ((TrF^(x)${x, x')F Xp (x')$ (x\ x)))) , (4) 



N c " " " ^ \ 2 Js Js 

where S is the minimal surface bounded by the contour C. Keeping only the string-generating 
field strength correlators, we have 

((TrF^(x)$(x,x')F Ap (x')$(x / ,x))) = (5, x 5 up - 5 pp 5 vX )D((x - x'f). (5) 
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In what follows we distinguish between the electric and magnetic contributions in Eq. (jSj), so 
that the structure da 0i d<Joi enters Eq. (J4]) multiplied by D , whereas the spatial part da^da^ 
is accompanied by D . The electric and magnetic string tensions are defined then according 
to Eq. flrj. 

We synchronise the quarks in the laboratory frame, putting x\q = x 2 q = t, and adopt the 
straight-line ansatz for the minimal string writing 

d<v(a;) = e ab d a w^(t, p)d b w u (t, p)dtdp, {a, b} = {t, P}, (6) 

where ^ t ^ t m3X , 0^/3^1, and the profile function is defined by the trajectories of the 
quarks, w^t, P) = fix^ + (1 — P)x 2 ^. For further convenience let us introduce two vectors: 

f=xi-x 2 , p = [(xi - x 2 ) x (/3xi + (1 - P)x 2 )] = ru, (7) 
which allow one to write the differentials in a compact form, 

daoi(x)d<Toi(x') = r(t)r(t')dtdt'dj3d/3', d<Tjk(x)d(jjk(x') = 2p(t) p(t')dtdt' d/3d/3' . (8) 
Presenting the averaged Wilson loop as 

-L(TrW(C)) = e- J , J = J E + J H , 

one can write for the electric and magnetic contributions separately: 

J E = dtdt' dp d/3' f(t)f(t')D E ((x-x') 2 ), (9) 
Jo Jo 

J H = / dt dt' / d(3 d(3' p(t, (3)p(t\ (3')D H ((x - x'f). (10) 



The correlation functions D E,H decrease in all directions of the Euclidean space-time with 
the correlation length T g which is measured on the lattice to be rather small, T g w 0.2 ^-0.3 fm 
[13J. Therefore, only close points x and x' are correlated, so that one can neglect the difference 
between r(t) and r(t'), p(t, (3) and p(t',j3') in Eqs. ffTU]) and also write: 

(x - x' f = (x(t, (3) - x(t', P')) 2 = g ab U b , i a = t-t\ i b = - P'. (11) 

The induced metric tensor is g ab = g a 5 ab , g 1 g 2 = det g = r 2 + p 2 = r 2 (l + u 2 ). Now, after an 
appropriate change of variables and introducing the string tensions, according to Eq. (jTJ), one 
readily finds: 

J E = a E I dt I dp = a E r / dt dp (12) 



y/r 2 + p 2 " Jo Jo VI + oj 



/■tmax rl „2 /"tmax rl , ,2 

J H = a H dt dp J. = a H r / dt / dp (13) 

Jo io Vr 2 + p 2 Jo Jo V 1 + ^ 2 

For <je = &h = <y the sum of J E and J H reproduces the well-known action of the Nambu- 
Goto string. 

For further analysis we resort, as was done in Ref. [12], to the Hamiltonian description of 
the quark-antiquark system under consideration. We also turn over to Minkowski space-time. 
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Then the Lagrangian of the quark-antiquark system can be derived from the exponent in 
Eq. ([3]) in the form: 



I 7£ / 7£ f 1 1 

L = —mi y 1 - - 777.2 Y 1 — x 2 — o^r I df3 

Jo 



yjl- [nx + (1 - /3)f 2 )] 2 



i 



+ a H r / dp = , (14) 

Jo ^1- [nx (^ + (1 - /5)f 2 )]2 r 

Two particular cases of the Lagrangian f[T4"j) are of most interest. The first such case corresponds 
to equal masses, whereas in the other case one mass is assumed infinitely large. Then, using 
the standard technique, one can proceed to the Hamiltonian of the quark-antiquark system, 

u lfp 2 r + m 2 \ [\Jcry v \ P 

H=- — hu + / dpi— 1 Yo 2 r H , , (15 

where £ = 1 for the case of equal masses (mi = m 2 = m) and £ = 2 for the case of the 
heavy-light system (mi — > oo, m 2 = m). Here <j\ = a H + rf^au — o\e), o"2 = ^]{pe — °u)- 
The fields //, v((3), and 77(/3) are the auxiliary fields, also called in the literature the einbeins 
[14^]. Generally speaking, the einbein fields appear in the Lagrangian and, even in absence 
of the corresponding velocities, they can be considered as extra degrees of freedom introduced 
to the system. The einbeins can be touched upon when proceeding from the Lagrangian of 
the system to its Hamiltonian and thus they mix with the ordinary particles coordinates and 
momenta. Besides, in order to preserve the number of physical degrees of freedom, constrains 
are to be imposed on the system and then the formalism of constrained systems quantisation 
[15J is operative (see, for example, Ref. [16] for the open straight-line QCD string quantisation 
using this formalism). A nontrivial algebra of constraints and the process of disentangling the 
physical degrees of freedom and non-physical ones make the problem very complicated. In the 
meantime, a simpler approach to einbeins exists which amounts to considering all (or some) 
of them as variational parameters and thus to taking extrema in the einbeins either in the 
Hamiltonian or in its spectrum [T7]. Being an approximate approach this technique appears 
accurate enough (see, for example, Ref. [18]) providing a simple but powerful and intuitive 
method of investigation. In this Letter we follow the given technique, so that extrema in all 
three einbeins are understood in the Hamiltonian (1151) . Notice that, for <je = &h — the field 
77 drops from the Hamiltonian and the standard expression for the string with quarks at the 
ends [12] readily comes out from Eq. ffl5|) . 

Notice that the kinetic part of the Hamiltonian ffl5|) has a very clear structure: the radial 
motion of the quarks happens with the effective mass /1, whereas for the orbital motion the 
mass is somewhat different, containing the contribution of the inertia of the string. 

We take the extrema in v{(3) and f](j3) now, approximating r){f3) by a uniform in f3 distri- 
bution. Then the Hamiltonian (Tl5|) takes the form 



2 The einbein [i is introduced in the Lagrangian via the substitution 2V AB — ► A/ [/,-{- Bfj, and allows one to 
simplify the quark kinetic term. The continuous einbein v(/3) enters through the same trick for the string term. 
Finally, the second continuous einbein rj(f3) appears due to the substitution A 2 / B — > —Brj 2 + 2r/A. As soon as 
extrema are taken in all einbeins, the initial form of the Lagrangian is restored. 
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and the spin-independent potential reads: 



Vsi(r) = i] a E r + ( — - rjo] °Hr + ^fiy 2 , r) = — -. — , (17) 
\?7o / q arcsm?/ 

with y being the solution of the transcendental equation 



ffl = l( 1+llS (I (18) 

a H r 2 Ay \ \ a H J J \r) J o H r 

The interested reader can find the details of a similar evaluation performed for the Hamiltonian 
(USD with a E = a H = a in Ref. [HI [20]. 

The remaining einbein fi is to be considered as the variational parameter to minimise the 
spectrum of the Hamiltonian ffTBT) . Obviously, the extremal value of \x depends on quantum 
numbers and acquires two contributions: one coming from the current quark mass m and the 
other, purely dynamical, contribution coming from the mean value of the radial component of 
the momentum p r . It is instructive to pinpoint the difference in the potential (1171) below and 
above the T c . 

At small r's, the potential (fTTj) turns to the centrifugal barrier /(/ + l)/(£/ir 2 ), whereas its 
large-r behaviour differs dramatically for the temperatures below and above the T c . Indeed, 
the leading large-r contribution to the inter-quark potential corresponds to ?/ < 1 and, for 
T < T c , reads: 

Konf(r) = a E r. (19) 

This is the linear confinement which is of a purely colour-electric nature and which admits 
angular-momentum-dependent corrections (see Refs. [T2[ [T9]). 

In the deconfinement phase, at T > T c , the colour-electric part of the potential ( 1TT1) van- 
ishes, the leading long-range term coming from the angular-momentum-dependent part of the 
interaction: 

K slW = f±iU... (20) 

Interestingly, in the deconfinement phase in absence of the confining potential, the spin- 
independent interaction becomes short-ranged decreasing as 1/r 3 at large inter-quark separa- 
tions. This feature means the full compensation of the centrifugal barrier which would naively 
behave as 1/r 2 instead. The reason is obvious: at large inter-quark separations, the effective 
quark mass fi is to be compared to the "mass" of the string or. The bound-state problem solved 
in the potential (|T9|) gives a large value (p r ) oc 0\e(?"), so that, even for light (massless) quarks, 
their effective mass fi appears quite large (/i ^> m). On the contrary, for light quarks and in 
absence of the strong confining interaction f|19p . the values of /i are small (/i m m) and can be 
neglected as compared to the string contribution cr^r. This makes the spin-dependent terms 
in the effective inter-quark interaction important in this regime, as opposed to the confine- 
ment phase, where they give only small corrections to the bound states formed in the confining 
potential f[T9~|) . 

In the next section we turn to the derivation of spin-dependent contributions to the inter- 
quark potential. 



3 Nonperturbative spin— dependent interactions 

In this section we return to the Green's function of the quark-ant iquark system ([3]) and restore 
spin-dependent terms. To this end we notice that the interaction of the quark spins with the 
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background gluonic field is to be added at the exponent. It enters in the standard combination 
VftuF^u, with cr^j, = ^(7^7^ —7^7^), and appears under the integral in the proper time r. After 
averaging over the background field, the resulting expression for the Wilson loop can be again 
written in the form of Eq. (J3J) , 

' (TrW(C)) = exp (-\ [ dn^{x) f dir Xp {x'){{TrF^{x)${x, x')F Xp {x')${x' , x))) V (21) 



but with the differential dn^ containing the spinor part, 

dn^ v (x) = ds^ u (x) - ia^dr. (22) 

The nonperturbative spin-dependent interaction appears from the combination of differen- 
tials involving a^ v . We skip the details of the derivation which can be found, for example, in 
Refs. [21| [22] and quote the result here. The leading spin-dependent term is the spin-orbit 
interaction (we omit contributions of Df ,H which bring about short-range terms 0(r~ 3 )), 

T/ (SiLi S 2 L 2 \f 1 dV 2 dV l \ 

Vso{r) = \-M-^)\rl^ + r^)^ (23) 
where V^(r) can be expressed through the colour-electric and colour-magnetic field correlators, 



^ = i r dT r ^rA), ?=-ir*f4i-^"M). (24) 

r dr r J J r dr r J J \ rj 



Notice that this result [2TJ [22], is not due to the 1/m expansion, but is obtained with the 
only approximation made being the Gaussian approximation for field correlators. Accuracy of 
this approximation was checked both at T = [23] and at T > T c [21] to be of the order of one 
percent. 

Only the V\ potential survives above the T c . Besides the corresponding /i-dependent de- 
nominator is to be corrected according to the discussion of the previous section — namely, one 
of /i's is to be augmented by the string rotation term yielding 



where 



V ^ = V(fr + 2W-f/2)W . dXD (T ' X) I 1 -r) ' <25) 

(v((i - ilif) = f dPv{P - m? = £fi (1 + >i) (~ - ■ (26) 

Jo z $ y Wo j 

and S = Si + S 2 , for the light-light system, and S is the light-quark spin, for the heavy-light 
quarkonium. 



4 Bound states of heavy quarks above the deconfinement 
temperature 

In the previous sections we derived the effective nonperturbative inter-quark potential, including 
spin-independent terms and the spin-orbital interaction. It follows from Eq. f[25l) that above the 
deconfinement temperature, for the states with the total momentum J = 1 + S, SL > and the 
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Figure 1: The profile of the effective potential (128]) for m = lGeV (solid line), m = 2GeV 
(dashed line), and m = 3GeV (dotted line). 



potential Vso{r) becomes attractive, with a possibility to maintain bound states. Furthermore, 
its slow decrease as r — > oo suggests that an infinite number of bound states exists, with the 
binding energies asymptotically approaching zero. Let us study these bound states in more 
detail. Hereafter <je = and we use the notation a for the magnetic tension an- 

In view of an obvious similarity of the light-light and heavy-light cases (the difference 
manifesting itself only in numerical coefficients), we investigate numerically only the light-light 
system, as a paradigmatic example. Furthermore, for r ^> T g , the potential fl^o^) does not 
depend on the form of the correlator D H since 



2 [ dr [ dXD H (r,X) (1 - - ) « o. 

J0 J0 V rJr>T g 



(27) 



Finally, we neglect the perturbative part of the inter-quark interaction for it is screened to a 
large extend contributing to short-ranged forces only whereas the effect discussed in this work 
is essentially a long-ranged effect. 

Therefore we study the spectrum of bound states in the potential 



/ arcsin y y \ 9 al 

V(r) = or + ay 2 — — j—-^, 28 

V y arcsin?/; w jxr(ju + 2(v(J3 - 1/2) 2 )) ' 



which is the sum of the spin-independent term ( fl71) and the spin-orbital term (125jl : y is the 
solution of Eq. ( fl8l) with oe = 0. In Fig. CD we plot the effective potential ( 1281) for three values 
of the quark mass: m = lGeV, 2GeV, and 3GeV. 

The resulting eigenenergy e nr i(fi) is added then to the free part of the Hamiltonian (TX5j) . 



2 

m 



M nr i (jj) = — +fi + e nA (ji) , (29) 



and this sum is minimised with respect to the einbein fi, 

= 0, M nrl = M nrl (/Jo). (30) 



dM n M 
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Parameter 


rrib, GeV 


m c , GeV 


m s , GeV 


a, GeV 2 


T g , fm 


Value 


4.8 


1.44 


0.22 


0.2 


0.2 



Table 1: The set of parameters used for the numerical evaluation. 





bb 


cc 


ss 


n r = 


-0.007 


-0.19 


-45 


n r = 1 


-5xl0" 4 


-0.015 


-2.7 



Table 2: The binding energy E Hr i = M nr i — 2m (in MeV) for the ground state and for the first 
radial excitation in the potential (I28I) with I = 1 for the bb, cc, and ss quarkonia. 



In Table 1 we present the set of parameters used in our numerical calculations, whereas in 
Table 2 we give the results for the binding energy for the bb, cc, and ss quarkonia above the 
T c for I = 1 and n r = 0,1. We ensure therefore that for I ^ the potential ( 1281) does support 
bound states. The binding energy is small (|-B nr /| <C T for the b and c quarks and \E Ur i\ < T 
for s quarks) so these bound states can dissociate easily. 



5 Bound states of light quarks above the deconfinement 
temperature 

In this section we turn to the problem of binding of light quarks. The effective potential (128)) 
admits different forms at different inter-quark separations, depending on which contribution, 
of the quark mass term // or of the "string mass" 2{v{(3 — 1/2) 2 )), gives the dominating con- 
tribution, that is for fi ^> or and \i <C or. If large distances contribute most to the bound 
state formation (the latter case), then V(r) = 0(1(1 + l)/(crr 3 )) + (9(//(/ir 2 )), where the first 
term comes from the spin-independent interaction (see Eq. (1241) ) and the other stems from the 
spin-orbit potential. The dependence of the binding energy on jj, is expected then to be rather 
moderate, approximately as 1//U. 

On the contrary, in the former case with the string dynamics giving a correction to the 
quark mass term, the potential fl28|) can be approximated as 

v(r) J±±A.^ (31) 

that is by the sum of the centrifugal barrier and the attractive Coulomb-like potential with the 
effective coupling 

al 

Ciett = —z- (32) 

The corresponding eigenenergy can be found in any textbook in Quantum Mechanics and gives 
a stronger dependence on \i, 

a 2 l 2 

e(n) oc -/i« 2 ff oc —. (33) 

Let us consider the states with / = 1 and n r = 0. We follow now the procedure described 
in detail in the previous section, solve the full problem numerically, and find the dependence 
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m [GeV] m [GeV] 

12 1 




Figure 2: The binding energy of the quark-antiquark system versus the mass of the quark for 
I = 1 and n r = (first plot) and n r — 1 (second plot). 



of the eigenenergy e on the einbein jj, to be 

£ M K -^79- ( 34 ) 

Comparing this to Eq. (1331) we find a good agreement, with the small deviation in the 
power resulting from the proper string dynamics. We conclude therefore that the dynamics 
of the system develops at the inter-quark separations T g <C r < m/a (since, for the set of 
parameters given in Table 1, m > aT g then there is room for such separations). 

Let us discuss now the procedure of minimisation of the spectrum (1291) in \i. First of all, 
let us notice that, in the einbein field formalism, the calculation of the spectrum naively looks 
like a nonrelativistic calculation due to the "nonrelativistic" form of the kinetic energy in the 
Hamiltonian with the einbein field /i introduced. In the meantime, the full relativistic form of 
the quark kinetic energy is readily restored as fi takes its extremal value and hence this is the 
procedure of taking extremum in /i in the masses (1291) to sum up an infinite series of relativistic 
corrections and thus to restore the relativistic spectrum. For example, the relativistic ground 
state eigenenergy Eq = m^Jl — (Za) 2 of the one-body Dirac equation with the Coulomb 
potential —Za/r can be reproduced exactly with the help of the einbein technique. Finally, 
one can visualise the form of /i considering the effective Dirac equation for the light quark in 
the field of the static antiquark source. When written in the form of a second-order differential 
equation, it contains the spin-orbit term of the same form as given in Eq. (T2"3]) but with \i 
replaced by the combination e + m + U — V, where U and V are scalar and vector potentials, 
respectively. For light (massless) quarks this combination takes drastically different values 
below and above the deconfinement temperature. Indeed, in the confining phase of QCD, when 
spontaneous chiral symmetry breaking leads to a strong effective, dynamically generated scalar 
potential U, this effective "/i" is large. On the contrary, above the T c , when U is small "//" is 
also small (it can be even negative since the eigenenergy e may have any sign). 

We find numerically that the extremum in \x for M nr ;(/i), Eq. (1291) . exists for m exceeding the 
value of approximately 0.22GeV (for the given a = 0.2GeV), and no extremum exists for smaller 
values of the quark mass (see Fig. [2] for the dependence of the binding energy E Ur i on the quark 
mass). This property of the bound state spectrum can be easily understood using the analogy 
with the bound state problem for the Dirac equation with the potential in the form of a deep 
square well or Coulomb potential discussed above. For example, for the Coulomb potential, a 



9 



problem appears as the coupling exceeds unity — the well-known problem of Z > 137. From 
Eq. (132]) we easily find this critical phenomenon to happen at m ss fi < y/a 0.4GeV. This 
estimate is in good agreement with the result of our direct numerical calculations quoted above. 

Physically this situation means that many quark-antiquark and/or gluon pairs are formed 
and finally stabilise the vacuum. Formally the problem is not anymore a two-body problem, 
but rather many-body, so that many-body techniques are to be applied. For example, in 
electrodynamics with Z > 137, one can derive the resulting self-consistent field of the Thomas- 
Fermi type [25J. A similar situation can be expected in the deconfinement phase of QCD. In 
absence of the linear potential, the einbein /j (playing the role of the effective quark mass) is not 
anymore bounded from below by the values of order y/aE — 0.4GeV coming from the binding 
energy in the linearly rising potential. To see the onset of this phenomenon in the framework of 
our two-body (one-body for the heavy-light case) Hamiltonian, one should take into account 
the negative-energy part of the spectrum, when the full matrix form of the Hamiltonian is 
considered [26J. Indeed, the matrix structure of the Hamiltonian occurs in the path-integral 
formalism from the two-fold time-forward/backward motion described by the positive/negative 
values of /x. Off-diagonal terms in the matrix Hamiltonian produce the turning points in the 
particle trajectory and result in Z-graphs. 

Notice that the same is true for the glueballs and gluelumps since in this case equations 
are the same as for light-light and heavy-light quarkonia, respectively, but with the quark spin 
replaced by the gluon spin and an by \oh- 

Concluding this section one can say that colour-magnetic (spin-dependent) interaction 
acting on light quark or gluonic systems enforces nonperturbative creation of light qq and 
gg pairs. 

6 Discussion 

The results obtained in this Letter allow us to comment on the general situation with the 
existence of bound quark-antiquark states in the deconfinement phase of QCD. First of all, 
contrary to naive expectations, the colour-Coulomb potential is screened down to a short- 
ranged interaction and bound states appear due to nonperturbative colour-electric (see Refs. [H 
[9]) and colour-magnetic interactions in the vacuum. Indeed, although quark-antiquark pairs 
in the relative S-wave cannot be bound by such interactions for T > 1.5T C (see, for example, 
Refs. [Bl El ED]), pairs with a nonzero relative angular momentum can form bound states at 
all T > T c since an grows with T. Second, formation of such bound states above the T c is 
energetically favourable since it lowers the system energy as compared to the ensemble of free, 
unbound particles. For heavy quarks the binding is weak and the system dissociates easily. 
Finally, the dependence of the binding energy on the quark mass is strong — the corresponding 
eigenenergy for strange quarks is around 10 -r lOOMeV rather then below lMeV for the charmed 
and bottom quarks (see Table 2). The situation becomes even more dramatic for the lightest 
quarks. The effective inter-quark potential for light quark flavours (and gluons!) becomes 
extremely strong and may lead to pair creation — the effect similar to the critical phenomenon 
in QED for the centre charge Z exceeding 137. Vacuum polarisation effects become important 
and they lead to a complete rebuilding of the vacuum structure of the theory. We anticipate a 
similar phenomenon to take place for light quarks (and gluons) in QCD in the deconfinement 
phase. 

It is important to notice that it was a separated quark-antiquark pair which was considered 
in this paper. In reality such quark-antiquark pairs are to be considered in the medium formed 
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by other quarks and gluons, that is as a part of the SQGP. As a measure of the interaction in 
SQGP one can consider the ratio of the mean potential energy to the mean kinetic energy of the 
particles in the plasma, T = (V)/(K). It is easy to estimate that (K) ~ T and (V) ~ cr H /T. 
This gives T = <jh/T 2 and so this parameter is large for quarks and it is several times larger 
for gluons. Therefore, SQGP is a strongly interacting medium which looks like a liquid, rather 
then as a gas. With the growth of the temperature the medium becomes more dense, and the 
mean distance between particles decreases. As this distance becomes comparable to the radius 
of the bound states discussed in this Letter, the latter will dissociate because of the screening 
effects. In other words, the hot medium plays the role of a natural cut-off for the effect of bound 
pair creation discussed above. Notice however that, for the quark masses around 0.2GeV, the 
radius of the bound state is of the order of one fm and it is expected to decrease further with 
the decrease of the quark mass, even if the pair creation process is properly taken into account. 
This means that indeed there is room for such bound states for the temperatures above the T c . 
Breakup, with the growth of the temperature, of such high-/ states for quarks, and especially 
for gluons which possess more degrees of freedom than quarks, may affect such characteristics 
of the plasma as its free energy and it entropy (for a recent attempt of explaining the near-T c 
behaviour of these characteristics see Ref. [27]). This work is in progress and will be reported 
elsewhere. 
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